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TROPOSPHERIC LIMITATIONS ON THE ACCURACY OF PHASE 
MEASUREMENT OF COORDINATES IN ASTRONOMY 


A. F. Dravskikh, A. A. Stotskiy, 

A. M. Finkel ' shteyn and P. A. Fridman 


This work considers the effect of tropospheric 
fluctuation effects on the accuracy of phase meas- 
urements of coordinates. The nature of the averag- 
ing of the trophospheric effects, if N coordinate 
measurements of duration T with period u are made, 
is Investigated. Various averaging modes depending 
on the relation of the various time parameters are 
investigated. Equations taking into account the 
correlations between individual observations are 
presented. It is shown that the correlation Inter- 
val between the individual observations is always 
greater than the fluctuation period of tropospheric 
inhomogeneities typical for a given baseline. 


In measuring the coordinates of celestial sources with the nelp /_ 0 
of ground-based optical and radio telescopes, the presence of tropo- 
spheric inhomogeneities of the refractive index leads to the genera- 
tion of random errors in the measured positions of the sources. If 
a double-element interferometer is considered as the telescope model, 
then the effect of the troposphere appears random advances in phase 
of the signal p in the troposphere along the routes to the first and 
second antennas, whose magnitude is determined by the fluctuation 

’'.Numbers in the margin indicate pagination in the original foreign 
text . 
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properties of the troposphere. As is well known, it is convenient to 
consider as a sufficiently general characteristic of these properties 
the structure function of the electric thickness of the tropo- 
sphere [1], which is related to the structure function of the phase 
in the geometric optics approximation in an elementary manner: 

/>*■• <!>(* + a/) - <i> (01*> *" (y) &!• 
where X is the wavelength. 


The time structure function of the real troposphere is well 
approximated by the Karman function: 






(i) 


with the parameters a =10 cm, t„ ."> h r and »; [2]. Here J'{.) is the 
gamma function and K, is the MacDonald function. The hypothesis of 
"frozen" turbulence, which is confirmed experimentally for the free 
atmosphere from the measurements of meteorological parameters up to 
scales the order of a kilometer [3], permits converting from the time 
to the spatial structure function />■(?.) where p is the inter- 

ferometer baseline and v is the average drift velocity of the tropo- 
spheric inhomogeneities , the average wind speed. The existing esti- 
mates for v are very indfinite [4]. However, if one assume a compro- 
mise value of v = 10 m/sec, which coincides with the typical wind /I 

speed for synoptic processes [5], then the overall scale of turbu- 
lence in the assumed model will be «•/» km. We note that 
inhomogeneities close to these scales have actually been observed [6]. 


In the region of scales small as compared to the overall scale 
of turbulence :••• (inertial interval), the Karman approximation (1) 
gives the well-known "5/3 law" following directly from the Kolmogorov- 
Obukhov turbulence theory [1]. For large values •.>•.„ (energy inter- 
val), the Karman function tends asymptotically to the constant value 
cu-ic-/)- 1 . determining the maximum advance of the phase difference in 
the troposphere. It should be noted that, although eke Karman 
approximation may not correspond to actual physical processes in the 
intermediate range of values of p [7], it remains a convenient 
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approximation, which is completely suitable for obtaining estimates 
In analytic form. 

The structure function permits evaluating the tropospheric 

error to arising with an instantaneous measurement of the coordinate 
of a source: to . This error can obviously be reduced if one 

uses- time accumulation. There is the widely held opinion that the 
nature of the averaging is determined only by the fluctuation period 
of tropospheric inhomogeneities pd* typical for the baseline p, 
which coincides with the time correlation interval of the fluctuations 
of the refractive index of the tropospheric inhomogeneities. It then 
seems that one can make .v T,t f independent readings during the time 
/ t f and, consequently, reduce the tropospheric error by \YT f times: 

i. '7~ 

so, » so . However, if one supports this point of view, then the fol- 
lowing paradoxical conclusion is Inevitable: the tropospheric angular 
error uo, in the inertial interval is smaller, the smaller the base- 
line. Actually, for p»- m .o p ” 1 (so—o, as — (l ) , and thus 
so, ^ v \\fi- » Fw’/v , since •*''/ — » as p — » . 

Actually, as rigorous analysis shows, in the general case 
does not determine the correlation interval and is always signifi- 
cantly greater than L forp<? 0 . The dependence of so, on the accumu- 
lation time T is determined by the mutual relation of the time /* char- 
acterizing the overall scale of the turbulence and the times >, and T 
characterizing the spatial and time fluctuations of the tropospheric 
inhomogeneities. In particular, in the inertial interval of greatest 
interest for optics, the time accumulation, as will be shown below, 
effectively weakens the effect of the troposphere only for suffi- 
ciently large accumulation times / ->s / n . Then iO, does not depend on 
the baseline. 

This work considers the nature of the time averaging of the 
tropospheric angular error and analytic estimates are obtained for 
the following model of observation procedure — N sequences of obser- 
vations of duration T with period ( : . y p- i) are carried o-t. The ob- 
tained results develop and generalize the results of [3] in which 
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some modes of averaging one observation of duration T were studied. 

In spite of the fact that a double-element interferometer is con- 
sidered as the model of the measuring device for analysis, the ob- 
tained results are also applicable to instruments with filled aper- 
tures. Actually, in the case of a continuous aperture with a uniform 
ample distribution, the fluctuations of the position of the principal 
lobe determined by the fluctuations of the center of gravity of the 
radiation pattern differ from the value obtained on the basis of 
fluctuations of the phase difference at the edges of the aperture 
only by a coefficient the order of unity (0.97 for a circular aper- 
ture [1], 1 - 1.5 for a linear aperture). The latter is related to 
the fact that small-scale fluctuations (the order of the dimensions 
of the aperture and smaller) are filtered by the continuous aperture 
and do not shift the center of gravity of the image in practice, 
but only lead to its blurring 

Spectrum of the phase difference fluctuations and general 
equations for the dispersion of the phase difference . For further 
analysis, it is convenient to transform from the structure function 
(1) to the phase fluctuation spectrum related by the Fourier 
transformation: 


which is related to the spectral density of the phase difference 
fluctuations Sy— *»- W at two points located at a distance p by the 
relation : 

**’i. (") *- '• 'in' jj- H'j, M. ( 3 ) 

where oj is the cyclic frequency. The periodic factor in (3) charac- 
terizes the spatial filtration of the tropospheric inhomogeneities 
with dimensions greater than the instrumental baseline p. 

Time averaging over the interval T leads to filtration of the 
fluctuations with frequencies <”>i IT : 

“Tv m *-* 
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The last expression for the spectral density determines the disper- 
b-.on Oj. the phase diiierence caused by motion of the tropospheric 
inhomogeneities for an instrument with a baseline p for a single 
sequence of observation of duration T: 


/% 
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where « - ? ? = T /2 , and 
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-et there be N sequences of observations repeated with period y 
w_„h a duration of the observations in each sequence . The 

time shift of the signal by the amount y is equivalent to multiply- 
ing -^ s spectrum by the quantity . Consequently, the total 
spectral density for N sequences will have the form: 
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Since the expression within the absolute value sign is a sum of 
terms of a geometric progression with the ratio , then, as can 
be shown: 

A -1 


I * -1 t ... + B ,v + 2^t l . ul |,v. 

-hus, the spectral density in the case of N measurements will be: 
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zjj integrating (3) by parts, we arrive at the following expression 
for the dispersion of the phase difference: 
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It is not difficult to show that the sum of the integrals 
*> = /» + /* = s »"- v . If A -W»-W»=° > then the individual sequences were 
averaged as statistically independent. Consequently, the group of 
terms related to the integrals A* ^ characterize the effect of 
the correlations between the individual sequences. 


and 


3y using simple trigonometric relations, the integrals /,. 
03 can be reduced to the sum of the Integrals of the form: 
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where S* («. l >. n ;) I s > 0 1 is the generalized hypergeometric series, while 
the integrals I 2 and can be reduced to the sum of the integrals: 
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which is the 3asset integral representation for the MacDonald func- 
tion [ 9 ]. 


Then, by keeping in mind the integral expressions (10), we have: 
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Here and below, (+) denotes summation over the corresponding 
two terms: / u • y) H* y)-t I U—u) . We draw attention to the obvious fact 
that as ; there is no correlation between the indi- 

vidual sequences and, as a consequence, a, =. a , \'X . 


The extremely cumbersome expression for the dispersion of the 
phase difference fluctuation (11) caused by the motion of tropo- 
spheric inhomogeneities becomes accessible for simple analysis in a 
majority of cases of greatest importance for observational astronomy. 

Asymptotic expressions for the dispersion of the phase dlffer - 
v :■ c . The nature of the averaging of a.*- depends on the relation of 
five time parameters: C ?'»-'-*• T~2ji and the total observation time 

= (.V-l) J * : T. . 

Cl 


A. Small t f . r e f Let all the characteristic time parameters 
t,. 7' e: -, ana consequently T, h , be small as compared to the time scale 
tQ determined by the overall scale of the turbulence. In other 
words, we will consider the case of small baselines ( K luo km) and 
small accumulation times ( r < i hr). 

Since the functions >’.(;) and A.., (*)(» >0) have the following 
asymptotic expressions for small * i : 
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we obtain for (11) after some transformation: 
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V.'e note that tne expression A -r^) determined by (13) 

coincides as : k — *■ and N = 1 with the analogous expression obtained 
by another method in [4]. We also draw attention to the fact that 
(A-r^U ,+ /.(r) = JV/ l (JW), as must be, since for *~T. r e f=‘ vr • 


In the approximation under consideration, the following rigorous 
inequality is satisfied for the terms characterizing the correlations 
between the individual sequences *?— 4* '.-j - r . . 


<»fSr3(- 4XA)”- 
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The right side of this inequality corresponds to the condition 

/ „</, when the effect of the correlation terms are minimum for 
e» 

given p, T and p. The left side of the inequality is satisfied 

C i 

for , when the effect of the correlations between the indi- 

vidual sequences are minimum for given p and T, and can serve as a 
good approximation of the exact expressions for 3 -> even for 
We also note that there is the estimate convenient for further 
calculations : 
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Figure 1. Dependence of 
>/ : j on the number of 
sequences N for T * 0.1 sec 

Dotted line — *r i» ; 
solid lines — - 1 '• * 

3 — 100 m 


Figure 2. Dependence of ** on 
number of observations N for 

; .r-. i". 

I . •- lou. ; |w, J - I *. r, Ml; T, l 

1 . •*' 

sec; dotted lines — 
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According to (14), the effect of the correlation terms are signi- 
ficant even for significant ratios :« /' and :»■*, (Figure 1). Actually, 
depends weakly on N and the simple accumulation of the number of 
sequences N for specified T e f and p does not lead to a significant 
cecrease in the tropospheric error (Figure 2). The latter actu- 
ally indicates that when all the characteristic time parameters are 
less than tQ, the effect of tine accumulation is scarcely perceptible. 


3. Small u and t p , large T e f. This mode differs from that con- 
side. i above in that the total observational time T e f can become 
larger than the characteristic time to for a sufficiently large num- 
ber of sequences N. In this case, as follows from (11), the summa- 
tion in expression (13) for 1^, J 2 and J 3 must be extended only to 
/.• A*-l, where N* is the Integer part of the ratio *•<> . Conse- 
quently, in this case, if , 
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An increase in the number of sequences in this mode can obviously 
^eaa to a significant decrease in 0 ; . 
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C. Small T and tp, large T e f.* Let the baseline p be small 
as compared to the overall scale of the turbulence pg and the accumu- 
lation time T in an individual sequence be less than the time to 
typical for the troposphere. If the Interval between Individual 
sequences ;* >/« , then <» , as is easily seen from (11), so that 
the sequences are statistically independent where, according to 
(12) and (13): 


- m . *i p«y 

Tv ^ i 0 -t- *)\ * / 


ikm 


(15) 


It is not difficult to see that within an individual sequence 
the time averaging occurs very slowly (T>i,) or not at all < O . 

This result is quite natural from the point of view of case A. 
Actually, the individual sequence of duration T can be considered 
as the sum of a large number of observations of duration T/N between 
which, according to (1*0, there are stroag ;orrelations. 


D. Small tp, large T. Let the accumulation time T in an indi- 
viuual sequence be greater than tg, and the baseline be small as 
compared to the overall scale of the turbulence. If ^ f is the 
time interval between the individual sequences, then we have beyond t 
dependence of the relation between AT, T and t p , according to (11): 
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(16a) 

(16b) 


When the time interval between individual sequences 
ana the individual sequences become statistically independent, which 
is quite natural since the characteristic time tQ determines the 


*Case C is the transition from case A to case B, since for ,.>*« 
1 and «, , -o , and for r e ^'*. . 
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natural scale of the correlation of the tropospheric inhomogenei- 
ties. Since , keeping in mind (12), we have according to (16a): 
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If the time interval sequences ±T<t t 
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Thus, for y>/„, the effect of the time averaging of the tropo- 
spheric advance of the phase difference becomes very perceptible, 
while the dependence of on the accumulation T is general and 

does not depend on the value of the index v. 


E. Large t p and small T. Let the dimension of the baseline p 
be larger than the overall scale of the turbulence pq, which occurs 
for observations on radio interferometers with extra long baselines 
(ELBR). The accumulation time for ELBR is limited by the instability 
of the time and frequency standards, and is always less than tQ. 


If- t 0 < in the approximation under consideration, 
according to (11) and (12): 


then 


/, ». 






and the individual sequences thus have maximum correlation since: 



(19) 


It is obvious that (19) is a direct consequence of (1) with ? > , 
and is the maximum possible advance of the phase difference in the 
troposphere. 


If for a small accumulation time T < L , the period between 
sequences p exceeds the characteristic scale to, then beyond the 
dependence on the relation between t P and p, the individual sequences 
become statistically i ;ependent, <>J o and 
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bosh cases, as should be expected, the tropospheric advance of 
the base difference depends neither on the duration of the sequence 
nor on tae baseline. 

Tropospheric limitations with coordinate measurements . The 
tropospneric advance of the phase difference »» t the expressions for 
which were obtained above, leads to a mean-square error in the meas- 
surement of the angular position of the source 1% = 3, > (radian) 

-0 1 **'•;, ii.i (arcsec) . 

The nature of the dependence of ao v ' on the baseline and the 
averaging time T e f is determined by the combined effect of two 
effects — spatial (the parameter t p ) and time (the parameters T, 
n ana T e p) filtrations. It is seen directly from the expression for 
the spectral density ( 8 ) that the effectiveness of both effects de- 
pends mainly- on the relation of V r > and the time parameter tg 
typical for the troposphere. The spatial and time averaging of 
- > becomes rather significant only when the corresponding time para- 
meter is greater than to* This means that the main contribution to 
the tropospheric error for measuring the angular position of a source 
in introduced by large-scale (low-frequency) tropospheric inhomo- 
jeneities. Their effect can be reduced only on baselines large as 
compared to the overall scale of the turbulence or for averaging 
times T large as compared to tg. 

1 . Baselines small as compared to the overall scale of the 
turbulence ?<*•(<,<*•). In this case, in accordance with the comments 
made above, a decrease in the tropospheric error can be achieved only 
a; to an increase in the accumulation time in the individual se- 
^t-nce or period tween sequences u, since the effect of the spatial 

filtration is sc ely perceptable in this mode. 

a) Accumulation time in the sequence T < /„ . If the period 
between sequences :* >d , then the observations in the individual 
so aencos are statistically independent and uo, = ao,, v A' > where, accord- 
ing to ( 13 ), the asymptotic expression for uo, will be: 


M, 


a::., >' *. 
1:1 r '\ 


If 

If 


rw, ' 

ryw,. ' 


( 20 ) 


where p is expressed in centimeters and T in seconds. The precise 
dependences of ao, on p and T are presented in Figure 3 and The 
dashed lines here indicate the results of calculations under the 
assumption that t, determines the correlation Interval between 
instantaneous measurements, while v TV is the number of independent 
readings during the time T. It is not difficult to see that this 



f, cm 


Figure 3* Dependence of the 
angular error of a single 
measurement a *i of duration T 
on the baseline S for 


VfM f* 



Figure Dependence of the 
angular error of a single meas- 
urement *’» over the baseline 
p on the accumulation time T 
for •, < 


point of view is true only in the trivial case T </,. , while for /lit 

T >* /, this leads to estimates differing by an order of magnitude 
from the actual. 


The mode T <i, , which corresponds to an instantaneous readout 
for a given baseline, is realized e; n emely rarely in practice (in 
radio astronomy — for large antennas and for bright sources) 
since the accumulation time T is determined not only by the fluc- 
tuation properties of the troposphere, but also by the noise proper- 
ties of the recording apparatus. As a rule, the signal/noise ratio 
becomes significant only for T — i, , and the tropospheric error (20) 
in this case does not depend on the baseline and decreases gradually 
with increasing T [8], 
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If the period between the sequences :* < l \ , then, as has 
already been stated, the individual sequences are not statistically 
independent, in connection with which the averaging of ao, in the 
general case goes more slowly than the law . According to (1*0 
and (i*4a), the asymptotic expressions for the angular tropospheric 
error have the form: 




if 

I* T > »,; 



if 



where p is measured in centimeters and u and T in seconds. 


( 21 ) 


Figure 5 presents the exact 
dependences of on the period be- 
tween sequences for the two modes 
(21). The estimates for other cases 
can be obtained by using Figures 
1-3. 

It is not difficult to see that 
the averaging law (21) coincides 
with the averaging law of statisti- 
cally independent sequences only 
for significant ratios •• /'• :• :• or 

v .v* . We note that for T v > 

ef 

and specified p, T, and u , one can always find such a value N, begin- 
ning with which an increase in the number of sequences does not lead 
to a notable decrease in the tropospheric error, and its dependence 
in this case on u has the general nature : AO * . 

Finally, the completely obvious consequences of (1*4) and (15) 
are the absence of a time averaging effect for T e f , » - M '~ ;r - : * •*' * . 

b) The accumulation time in the sequence T>U . in the mode 
when the accumulation time T is greater than the characteristic tQ, 
there occurs time averaging of large-scale tropospheric inhomogeneities 



Figure 5. Dependence of 
the angular error a** on 
the period between sequences 
for T - 1 sec: 

1 — N = 2; 2 — 5; 3— 

10 ; *4 — 50 ; 3 — 100 
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(with frequencies - ,• i ) which, as has already been noted, make 
the .tain contribution to the error in measuring the angular position 
Of a source. If the time interval between sequences sr~?-T < /. , 
then according to ( 18 ) the averaging of the tropospheric error will 

be very intensive: 

. |H* - - 

— rr ~\ i (22) 

where the more effective, the smaller AT (Figure 6). The latter is 
quite obvious since the case of continuous measurements o ( T ST) 

6a 

is the most favorable for decreasing the tropospheric error. 

If then individual sequences become statistically inde- 

pendent, and according to (19): 


ik, r* • i'** * V r. 


(23) 


The tropospheric error does not depend on the baseline p in both 
cases under consideration. 


We note in conclusion that 
even for a very large total obser- 
vation time, according to (22) and 
(23), it is difficult to obtain 
accuracy better than 0.05 arcsec 
(Figures 3, 4, 6). 

2. Baselines large as com- 
pared to the overall scale of 



Figure 6. Dependence of 
the error of the angular 
measurement on the time 
interval between sequences sr 
for T * 10 sec 


the turbulence f- > P. <<, > *•) . In 
this case, which is typical for 
measurements of ELBR, T<t, , and there is no time averaging within 
the sequence since: 


(24) 


where p is measured in centimeters. A decrease in the tropospheric 
error occurs only due to saturation of V in the energy interval, 
and the effect which in some sense "clears" the troposphere. In 



spite of th 


m. 


.HUITY QV 


b wo • f w °&5? 

nRlC.FiA' ' ' 


JIODUC’ 

ft$e\ qCi^MritfCwe raging, the effect of tropospheric 


effects can be limited with the help of ELBR on baselines the order 
of several thousand kilometers, as is well known, to thousandths of 
arcsec per unit measurement. According to (2), ELBR observations 
can be considered statistically independent only for i» when 


A A \ ,v . 


Conclusion . The presented estimates show that the fluctuation 
effects of the troposphere strongly limit the accuracy of coordinate 
measurements even for large accumulation times (7* > <•) . This conclu- 
sion is true to the greatest degree for measurements carried out on 
small baselines (?<?.) , which are typical for optical astronomy. 

The troposphere limits the accuracy of the coordinate measurements to 
hundredths of an arcsec even with exposures of several hours. 

The use of the methods of differential astronomy or the refer- 
ence object method can significantly weaken the effect- of the tropo- 
sphere. However, this assertion is true only for very small angular 
distances $ between the observed and reference objects: p „=$h ^,<p, 

where h ei * £ 10 km is the effective height of the troposphere; p is 
the diameter of the telosoope or the baseline of the Interferometer. 
In this case, the phase advance is determined by the effective base- 
line P e: -, and the angular error will thus decrease with decreasing $. 
In connection with this remark, all the estimates obtained in the 
work are completely applicable both for absolute as well as for dif- 
ferential measurements in optical and radio astronomy. 
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